Ground-state properties of the one-dimensional attractive Hubbard model with 

confinement: a comparative study 



O 

(N 

in 



C/3 



o 



o 



X 



Ji-Hong Hu,^ Jing-Jing Wang/ Gao Xianlong/'0 Masahiko Okumura,^'^ 
Ryo Igarashi,^''^ Susumu Yamada,^'^ and Masahiko Machida^'"^ 

''Department of Physics, Zhejiang Normal University, Jinhua 340012, China 
^CCSE, Japan Atomic Energy Agency, 6-9-3 Higashi-Ueno, Taito-ku, Tokyo 110-0015, Japan 
''CREST (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan 
(Dated: January 6, 2011) 

We revisit the one-dimensional attractive Hubbard model by using the Bethe-ansatz based density- 
functional theory and density-matrix renormalization method. The ground-state properties of this 
model are discussed in details for different fillings and different confining conditions in weak-to- 
intermediate coupling regime. We investigate the ground-state energy, energy gap, and pair-binding 
energy and compare them with those calculated from the canonical Bardeen-Cooper-Schrieffer ap- 
proximation. We find that the Bethe-ansatz based density-functional theory is computationally 
easy and yields an accurate description of the ground-state properties for weak-to-intermediate in- 
teraction strength, different fillings, and confinements. In order to characterize the quantum phase 
transition in the presence of a harmonic confinement, we calculate the thermodynamic stiffness, 
the density-functional fidelity, and fidelity susceptibility, respectively. It is shown that with the 
increase of the number of particles or attractive interaction strength, the system can be driven 
from the Luther-Emery-type phase to the composite phase of Luther-Emery-like in the wings and 
insulating-like in the center. 

PACS numbers: 71.15.Mb,03.75.Ss,03.75.Lm,71.10.Pm 



I. INTRODUCTION 

The Hubbard Hamiltonian is a simplified but im- 
portant prototype model for strongly correlated elec- 
trons in solid state. The one-dimensional (ID) version 
is usually used as a basic model examining the elec- 
tronic correlation behavior in the quasi-one-dimensional 
systems exhibiting the itinerant character of electrons 
and superconductivityii^ like the organic superconductor 
(TMTSF)2X,i^ Due to its solvability by the Bethe-ansatz 
method^ this model serves as a popular benchmark to 
test the accuracy of different approximate methods in 
handling the ground-state (GS) properties of strongly 
correlated fermions especially in an intermediate cou- 
pling strength!^""— This model has recently arisen renewed 
interest, in ultracold atomic gases, where fundamental 
many-body phenomena can be investigated in a well- 
controlled way. 

The ultracold atomic gases loaded on optical lattices 
have opened an exciting field in simulating condensed- 
matter model Hamiltonians like the boson and fermion 
Hubbard models. More recently a Tonks-Girardeau gas 
of bosonic ^^Rb atoms was realized experimentally in a 
ID optical lattice)^ In the case of fermionic atom gases, 
however, it is more difficult to cool them down because 
of the Pauli exclusion principle different from the bosonic 
atoms. Up to now, a two-component Fermi gas of ^'^K 
atoms with tunable interacting strengths has been pre- 
pared in a quasi- ID geometry;^ An interacting Fermi gas 
of ^"^K atoms has been demonstrated in three-dimensional 
optical lattices.— Within present-day techniques,— it is 
possible to trap interacting gases of fermionic atoms 
in ID optical lattices and then simulate the ID Hub- 



bard model experimentally. At the same time, experi- 
mental achievements on population-imbalanced ultracold 
fermion o^^d^ make it possible to systematically study the 
exotic paring states in ID fermion systems with or with- 
out optical latticesJ^ Moreover, the square- well trap and 
optical box-like trap have been experimentally realized in 
the atomic system, which makes it meaningful to study 
the physics of ID uniform atomic gases with a hard wall 
directlyji^ 

Considering the strongly correlated nature of ID sys- 
tems and the enhanced quantum fluctuations, it is crucial 
to have techniques beyond the mean field theory, both 
analytically and numerically. Among them, Bethe-ansatz 
solution^ quantum Monte-Carlo (QMC) simulation^"— 
density-matrix renormalization group (DMRG)r^"— and 
some sophisticated theoretical methods, such as the dy- 
namical mean field theory combined with the numerical 
renormalization group^^ and a time-evolving block deci- 
mation numerical study^i, are listed as well-established 
techniques. On the other hand, the exact diagonaliza- 
tion scheme^^i^ has also been used to obtain the exact 
numerical solution for the small size systems as an ulti- 
mate technique. 

The inhomogeneity caused by the boundary condi- 
tions or the external potentials, normally invalidates a 
reliable analytical method that is usually used in the 
homogeneous system. On the other hand, the numer- 
ical schemes mentioned above, e.g., the exact Bethe- 
ansatz solution incorporated with local density approxi- 
mation is proved to be useful in obtaining the phase di- 
agram and the collective oscillations of the atomic mass 
density in the presence of harmonic traps The lat- 
tice version of the Bethe-ansatz based density-functional 
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theory (BADFT)^i22i"i2^ is also a good candidate in 
studying inhomogeneous strongly-correlated system. For 
the strongly interacting regime, some non-perturbative 
methods, like, the Fermi-Bose mapping and the Bethe- 
ansatz method, have to be usedi^i^"— In general, the 
Bethe-ansatz solution provides reliable results beyond 
the mean-field theoryi^ BADFT has been used in study- 
ing the ID Gaudin-Yang gases^^^ and ID Hubbard 
model)^'^ which is different from the density-matrix 
functional theory with the basic variable to be the single- 
particle density matrix4S 

The rich non-uniform phases induced by the trapping 
potential can be identified by the variance of the lo- 
cal density, the local compressibility)^ the double oc- 
cupancy;^ or the thermodynamic stifFness4^ The entan- 
glement entropy is found by Franga et al. to be a pow- 
erful measure characterizing the phase-separated states 
in spatially inhomogeneous environment.— Recently, a 
new concept named fidelity in quantum-information and 
quantum-computation theory, is used to characterize 
quantum phase transitions by computing the overlap be- 
tween two states of nearby points in parameter space. 
By measuring the similarity between these two states, 
it is expected to show a minimum in the fidelity due 
to a dramatic change in the ground-state wavefunctions 
around the critical point of the quantum phase transition. 
The studies are made both experimentally and theoret- 
ically on many model systems, such as on the ID XY 
model, the Dicke model, and the ID Hubbard model4^ 
Gu further proposed a density-functional fidelity to mea- 
sure the similarity between density distributions of two 
ground states in parameter space based on the obser- 
vation that the density distributions^ can reflect the 
change in the ground-state wave functions. An inter- 
esting question then arises: Can this method be applied 
to quantify the composite quantum phases transitions of 
the ID Hubbard model in the presence of confinements? 

The ID Fermi-Hubbard model with attractive on-site 
interactions, has long served a "minimal" model that best 
describes superconductivity. For the fermionic atomic 
gases loaded on optical lattices, the strength or even the 
sign of the on-site interaction between different species 
can be easily tuned by using a technique named Fcsli- 
bach resonance4^ Marsiglio calculated the ground-state 
energy (GSE) and energy gap to the first excited state 
for the attractive Hubbard model using the variational 
canonical Bardeen-Coopcr-Schricffer (CBCS) and grand 
CBCS wave functional This model for different con- 
centrations of electrons is also studied within a self- 
consistent field methodji 

In this work, motivated by earlier theoretical work and 
by ongoing experimental efforts, we revisit the ID at- 
tractive Hubbard model in the weak-to- intermediate cou- 
pling strength under different boundary conditions with 
or without trapping potential. We present a fully micro- 
scopic theoretical study on the ground-state properties. 

We begin with introducing the Bethe-ansatz solution 
of the model and solving the coupled equations. In Sec- 



tion mil we provide a Bethe-ansatz based spin-density- 
functional theory (SDFT) and its local-spin-density ap- 
proximation. In Section IIV| we show our main numer- 
ical results, compared to the exact Bethe-ansatz, the 
CBCS^, and DMRG ones. Our comparisons are done 
on general band fillings and a wide range of coupling 
strength in three different systems, i.e., ID lattices with 
periodic boundary condition, hard-wall, and harmonic 
confinement. The energy, the binding gap, and the pair- 
binding energy of the ground state are in a good agree- 
ment with the exact results. We also provide a simple and 
direct diagnostic of the quantum phase transition by cal- 
culating the thermodynamic stiffness, density-functional 
fidelity, and fidelity susceptibility. At last, our results are 
summarized. 



II. THE ID ATTRACTIVE FERMI-HUBBARD 
MODEL 

We consider a ID attractively interacting two- 
component Fermi gas with iVf atoms in a lattice with 
unit lattice constant and Ns lattice sites, which can be 
described by a one-band Fermi-Hubbard model,— 

Hs = Hycf + Hcxt 

i— 1,(7 i—1 
+ ^ ViaflicT , (1) 

Z— 1,(T 

where a =t,i is a pseudospin-1/2 label for two internal 

hyperfine states, fii = J^a^i'^ = J^a'^la^ia is the to- 
tal site occupation operator, t is the tunneling between 
nearest neighbors, U is the strength of the on-site at- 
tractive interaction, and Via- describes the trapping po- 
tential. We consider in this paper ID lattices with pe- 
riodic boundary condition, in a simple box-shaped trap- 
ping potential (Vicr = inside the box and Vi^ = +oo 
elsewhere, namely, a hard-wall) and a parabolic confine- 
ment {Via = V2[i- {Ns - l)/2]2 with V2 the amplitude 
of the confining potential) . The numbers of atoms with 
spin up and spin down are N-^ and N_i, respectively. 

The ID attractive Fermi-Hubbard model in the ab- 
sence of confining potentials (i.e., Hs = iJrof) belongs 
to the universality class of Luther-Emery liquids, char- 
acterized by a massive spin gap and zero charge gap. 
At zero temperature, the properties of Href in the ther- 
modynamic limit {Na,Ns — >■ 00 with finite N^/Ns) are 
determined by the filling = N^/Ng and by the dimen- 
sionless coupling constant u ~ U/t. According to Lieb 
and Wu^, the ground state of H^cf for different fillings 
in the thermodynamic limit is described by the coupled 
integral equations for the momentum distributions p{k) 
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(3) 

(4) 
(5) 



{u/2Y + (A - A')2 ■ 

The non-negative parameters Q and B are determined 
by the normahzation conditions d/c p(fc) = 1 — 2s 

and J'_^^ c?Acr(A) = (n — 2s)/2, where the magnetization 
s is defined as s = {N^ — Ni) /{2Ns). The ground-state 
energy per site of the system is written in terms of the 
momentum distributions as, 



eGs{n,s,u) = -\U\ 





'n 




\J\ 













dkp{k) cos k, (6) 



which is a function of filhng n, magnetization s, and inter- 
action strength u. The GS properties of the ID attractive 
Hubbard model have been analyzed by many authorsj^i^ 



III. LATTICE SPIN-DENSITY-FUNCTIONAL 
THEORY 

Density functional theory (DFT)^i^ is a powerful 
tool to calculate the GS properties of a many-body 
Hamiltonian including inhomogeneity. Through Kohn- 
Sham mapping, the true system of interacting particles 
is mapped onto a fictitious one in which the particles 
do not interact but reproduce the true particle density. 
The central problem in DFT is the exchange-correlation 
(xc) functional, which remains elusive and can be sim- 
ply treated by the so-called local-density approximation 
(LDA)j^ The essence of LDA is to locally approximate 
the xc energy of the interacting inhomogeneous system 
based on that of an interacting homogeneous reference 
system. In discrete lattice systems, there is a lattice- 
version DFT, the so-called site-occupation functional 
theory (SOFT), introduced in the pioneering papers by 
Gunnarsson and Schonhammer to study the band-gap 
problem in the context of ab initio theories of funda- 
mental energy gaps in semiconductors and insulatorsj^ 
Within SOFT the GS site occupation can be obtained 
by solving sclf-consistcntly the lattice Kohn-Sham (KS) 
equations 
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FIG. 1: (Color online) Top panel: The exchange-correlation 
potential (in units of t) of the ID attractive Hubbard model 
as a function of n for various values of s with u = —1. The 
upper three data are for V^°^ and the lower three for V;''^°f • 
Obviously there is no charge gap for the attractive Hubbard 
model. Bottom panel: The same as that in the top panel but 
for u = —4. 



together with the closure relation 



(8) 



Here the sum a runs over all the occupied orbitals, and 
the effective KS potential is given by T^j^^[7iiCT] = Un.i^ + 
V^'^l'Tiia] + V°^^ where a = -a. The first term in the 
Kohn-Sham potential comes from the Hartree mean-field 
contribution, while ^^^'^[j^io-] is the derivative of the xc 
energy i?xc[%, '^j.] evaluated at the GS site occupation. 
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FIG. 2: (Color online) Top panel: The exchange-correlation 
potential (in units of t) V^c"™ of the attractive Hubbard model 
as a function of s for various values of n with u — —4. The 
spin channel opens a gap at s = 0. Bottom panel: The same 
as that in the top panel but for V^^"™. 



FIG. 3: (Color online) The xc gap Axc,ct('t., it) (in units of t) 
as a function of |it|. In the positive direction, A^c,T-{n,u) are 
shown and in the negative direction Axc,4.(7^, u) are plotted. In 
the inset, the semilog plot for Axc,t("-5 ^) is shown for \u\ —0.1 
to 0.5. 



where the xc potential of the ID homogeneous Hubbard 
model is defined hy^ 

d 

V^°^{n, s, u) = - — [eGs(n, s, u) - ecsin, s, 0) - Un,„n,a] 

Here, tQsin, s,u) is the GS energy per site of the ID 
homogeneous system defined in Eq. ^ and e^^^in, s, 0) = 
— (4i/7r) sin(7rn/2) cos(7rs). 
Making use of the chain rule, 



d 

dUa 



On d ds d 



dua dn dua ds 



we obtain an explicit expression for the exchange- 
correlation potential in terms of the chemical potential 
/i(n, s, u) and magnetic field h{n, s, u), 



The total GSE of the system is given by 
Eo [nt, = E E e("' -EE 

- E E ^ ^i'y'^is + -^xc [n^, »i] ■ (9) 



The local-spin-density approximation (LSDA) for the 
spin-polarized system has been shown to provide an ex- 
cellent account of the GS properties over a large variety of 
inhomogencous systems. In this work we employ the fol- 
lowing Bethc-ansatz based LSDA functional (BALSDA), 



ia I BALSDA 



(10) 



V^::^{n,s,u) = /i(n,s,u)±-/i(n,s,u) 



-I- 2t cos 



f/ns, (11) 



where the upper (lower) sign refers to the spin-up (spin- 
down) atomic species. The chemical potential and mag- 
netic field are defined by the first derivative of the GSE 
of the homogeneous system with respect to filling n and 
magnetization s as. 



h{n, s, u) 



ds 



ecsin, s,u) . 



(12) 
(13) 
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is always a discontinuity at s = 0, which corresponds 
to a gap opened in the spin channel signalling a Luther- 
Emery liquid phase of paired particles. We define the xc 
gap Axc,cr as the discontinuity in V^°^{n, s,u) at s = 0, 

A,e,.(n,u)- lim V^°^{n,s,u)- lim V^°^{n,s,u). 

The results are shown in Fig. [3] for different fillings. 
We notice there is a symmetry between Axc,t(fi, u) and 
^xc,ii'n,u), as a result of the symmetries of the GSE 
^Gsin, s, u) = ecsin-, —s, u) and our choice for the Hartree 
energy. As a consequence, the xc potential in Eq. ((TU]) 
possesses a discontinuity at s = and naturally the infor- 
mation of the Luther-Emery nature and the spin energy 
gap in the reference system is transferred to the inhomo- 
geneous system through Eq. ((TU)) . 
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FIG. 4: (Color online) Top panel: The ground-state energy 
per site (in units of t) of the attractive ID Hubbard model 
with periodic boundary condition as a function of particle 
density n = N{/Na for « = — 2 in a lattice with A^s = 32 
sites. SDFT/BALSDA results (crosses) are compared with 
the CBCS data (triangles) and the exact results (open circles). 
The CBCS data are taken from Ref.|4^ and the exact data are 
calculated from the couple Bethe-ansatz equations. Bottom 
panel: The same as that in the top panel but for u = —4 in 
a lattice with As = 16 sites. The system of even number of 
particles is implicit with zero magnetization (s = 0). For the 
system of odd number of particles, we choose A-|- — Aj_ = 1, 
that is, s — l/(2As). In this case, SDFT has to be used to 
obtain an accurate GS energy. 



In Fig. [U we illustrate the numerical results for the xc 
potential as a function of filling n for different magneti- 
zation s. Obviously, the xc potential is continuous in the 
whole range of n and there is no charge gap for the sys- 
tem of attractive interactions, which is different from the 
repulsive case where a discontinuity occurs and a charge 
gap opens at n = li^ 

In Fig. [21 we plot the xc potential as a function of s 
for various values of n with u = —A. We see that there 



IV. NUMERICAL RESULTS 

In this section we present the numerical results ob- 
tained from the self-consistent solution of Eqs. ([7])- 
(0) based on the BALSDA of Eq. We illustrate 

our main numerical results, which are summarized in 
Figs. Hl lTUl We compare the results obtained with the 
SDFT/BALSDA to results of DMRG and CBCS. In the 
use of DMRG, the number of states kept is set maximum 
1200 and the truncation error is less than 10~^". 

First, in Figs. HHHl we show the results for systems of 
periodic boundary condition. In Fig. 31 the GSE per site 
(absolute value) is plotted as a function of the particle 
density for u = —2 and u = —4. Good agreement is seen 
in the whole range of fillings between the SDFT results 
and the exact Bethe-ansatz ones, while CBCS calcula- 
tions deviate in higher density. In Fig. [5l we show the 
GSE as a function of interaction strength for the system 
of even and odd number of particles. The comparison 
shows that SDFT results give poorer agreement for small 
system size with even number of particles than CBCS but 
better agreement in other situations, such as for larger 
system size or stronger coupling strength. We should 
keep in mind that the reasons for which the SDFT re- 
sults deviate from the exact Bethe-ansatz ones are due 
to both finite size effects and the xc functional built in 
the KS potential, while for the CBCS results they are 
affected by both finite-size effects and lack of correlation. 
The ground-state energy and the energy gap (in Fig. [51 
see below) are less accurate by SDFT in small size sys- 
tems and by CBCS at intermediate coupling strengths. 

For the ID Hubbard model of attractive interaction, 
one defines the gap or binding energy as^^i^ 

An, = \ [EoiNi - 1) - 2EoiNt) + EoiNf + 1)] , (14) 

which describes the energy difference between two sys- 
tems, the one of two subsystems with Nf particles each, 
and the other of A'f + 1 and A'^f — 1. respectively. From 
the viewpoint of DFT, this gap has two contributions. 
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FIG. 5: (Color online) The ground-state energy per site (in units of t) of the attractive ID Hubbard model with periodic 
boundary condition as a function of the coupling strength u with even (left panel) and odd (right panel) number of particles 
in a lattice of A'^s = 8 (upper frame) and = 64 (lower frame) sites for various values of fillings. SDFT/BALSDA results (x) 
are compared with the exact (solid line) and CBCS (dotted lines) ones, (a) For Ns — 8 and A'^s = 64 with even numbers of 
particles (the magnetization s = 0). (b) The same as (a), but with odd numbers of particles (the magnetization s — l/{2Ns), 
and Nt-Ni^ 1). 



the one from the KS gap, and the other from the xc gap, 
related to the discontinuity of GSE as the number of 
fermions varies across an integer <^ This definition is also 
explained as the negative of the excitation gap, which 
equals to the half of the energy it takes to break a pair in 
the system<^ In Fig. [51 the energy gap as a function of 
the density is shown for (a) u = —1, (b) u = —1.5, and 
(c) u = -A for TVs = 4 (left panel), and = 16 (right 
panel). The comparison shows the SDFT calculations 
are qualitatively right in describing the properties of the 
gap. Quantitatively we find that a better agreement is 



achieved by SDFT for larger system size and stronger at- 
tractive interaction. The gap is negative for odd number 
of particles, and positive for even number of particles, 
which reflects the fact that the system with negative in- 
teractions favors pairing4ii^i^ This is the even-odd ef- 
fects discussed by Tanaka and Marsiglio42^ From the 
right panel of Fig. [51 we notice that the energy gap oscil- 
lates as a function of even particle numbers (depending 
on whether the number of particles is 4m or 4m -I- 2), 
which is the super-even effect (for details, see Ref. lish . 

The confined system composed of the hard wall, which 
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FIG. 6: (Color online) The energy gap (in units of t) of the attractive ID Hubbard model with periodic boundary condition 
as a function of the particle density n = Nf/Ns is plotted for various coupling strength in a lattice of sites A^'s = 4 and 16. 
The SDFT/BALSDA (crosses) results are compared with the exact (open circles) and CBCS data (upper triangle). The left 
panel is for Ns = 4: and the right panel A'a = 16. The gap is negative for odd A^f but positive for even A'f reflecting the pairing 
information induced by the negative interactions, (a) For u — —1; (b) For u = —1.5; and (c) For u = —4. For the system 
of even number of particles, the zero magnetization (s = 0) is used. For the system of odd number of particles, we choose 
m-Ni = 1, that is, s = l/(2N,). 
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FIG. 7: (Color online) (a) The ground-state energy per site 
(in units of t) of the attractive ID Hubbard model in a hard- 
wall as a function of particle density n — N[/Ns in a lattice 
of A'^s = 32 sites for u = — 1 and —4. 



can be realized in experiments with the atom chips, fo- 
cused laser beams and trapped impurity atoms breaks 
the inhomogeneity but still can be solved through the 
Bcthe-ansatz method<^ In Figs. [7] and |S1 we show the 
SDFT results for the GSE compared to DMRG ones. The 
system is of Ns = 32 lattice sites with u = — 1 and —4. 
Wc find that, in almost all the range of particle density 
and wcak-to-intermcdiatc interaction strength, SDFT re- 
sults give an excellent agreement with DMRG calcula- 
tions. Small deviations happen around the half filling. 

Next, we discuss a ID system of Fermi gases trapped 
in a harmonic potential, which is a system that breaks 
the homogeneity and the integrability of the model and 
can not be exactly solved. A quasi-classical approach 
called the local density approximation is usually used to 
study the density distribution profile when the external 
potential is slowly varyingi^ Here and in the following 
the DFT based on LDA is used. We would like to mention 
that for the system of even number of particles and spin- 
unpolarized, the DFT is used. But for the system of 
odd number of particles, the SDFT has to be used. Our 
results are shown in Table HI and HIl and Figs. 1^ and [TUl In 
TableHl we show the GSE for small systems of lattice size 
Ns = 100 and large systems of Ng = 300 for the weak- 
to-intermediate interaction strength. The corresponding 
pair-binding energies are calculated accordingly, which 
are defined as the energy difference between two systems, 
the one with two subsystems of particles Nf and A^f -I- 2, 
the other with two subsystems of A^f particles each, 

£p{Nf) = Eo{Nt + 2) + Eo{Nf) - 2Eo{Nf + 1) . (15) 

Our calculation shows that the pair-binding energies arc 
negative, which illustrates that the system likes to accom- 
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FIG. 8: (Color online) The ground-state energy per site (in 
units of t) of the attractive ID Hubbard model in a hard-wall 
as a function of the coupling strength it in a lattice of sites 
A^s ~ 32 for various values of fillings n. 



modatc pairs and has the tendency to form stable dimers. 
Furthermore, we calculate the four-particle binding ener- 
gies;^ 



(16) 



where AEi{N{) is the energy difference defined as 
AE,{N{) = EoiNf + i) ~ Eo{Nf). A bound state of 
two particles is confirmed by observing that, in Tabic HIl 
£p{Nf) < and £i{Ni) > 0, which shows the evidence 
of a paired superfiuid. We conclude that, the system in- 
duces effective attractions between particles leading to 
pairs formation {£p{Nf) < 0) and effective repulsions 
between pairs (£^4(A^f) > 0) avoiding clusters forma- 
tion. The same calculation is done for the ID Fcrmi- 
Hubbard model with repulsive interactions. It is found 
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FIG. 9: (Color online) The thermodynamic stifTness Sp (in 
units of t) as a function of A't for it = — 1 and u = —2 with 
V2A = 10^'^, and TVs = 200 lattice sites. In the inset we 
show the density profile in the bottom (top) for Ni — 40 
(A^f = 132) for u = —2. Their corresponding positions are 
indicated by the arrows in the A^/-axis. The vertical line 
indicates another possible quantum phase transition driven by 
the attractive interaction strength at fixed A^f — 128, which 
is checked through the concept of density-functional fidelity 
and fidelity susceptibility in Fig. (TO) 



that £p{N{) > and £4{Nf) > and no paired superfluid 
is formed. The comparison of the GSE between SDFT 
and DMRG shows that the agreement is up to the second 
decimal digit, except for few cases, such as for N{ = 30 
and u = —2 with the relative error of 1.2%. 

The harmonically trapped system described by Eq. ([T]) 
allows for two different spatial coexistent quantum 
phases. The one is the phase characterized by the Luther- 
Emery (LE) liquid combining with the appearance of the 
atomic density wave (n^ < 2), and another is charac- 
terized by a composite phase with a plateau in the trap 
center of tightly bound spin-singlet dimers (rii = 2) sur- 
rounded by Luther-Emery layers (n^ < 2). We call the 
latter the band-insulating (BI) phase. These two differ- 
ent phases can be well classified by a scaled dimen sionless 
variable, the characteristic density p = Ni^Yijt and the 
interaction strength u in the p-u phase diagram^. Thus 
for certain u, the quantum phases in the presence of con- 
finement can be driven from one to another by increasing 
p. Here we study this phenomena by investigating the 
thermodynamic stiffness, which is defined as the inverse 
of the global compressibility. 



Sfj, 



5Nf 



where p is the global chemical potential. The stiffness has 
been successfully used in the confined system to char- 
acterize the quantum phase transitions. The different 
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FIG. 10: (Color online) Density function fidelity (upper 
panel) and fidelity susceptibility (lower panel) as a function of 
couple strength u for = 100, Ni = 128, and ¥2/1 = 10~^. 
5u = 0.004 (solid red line) and Su = 0.002 (dotted blue fine) 
are shown in the figures. 



coexisted phases appeared in the harmonically confined 
system are well separated by the nonanalyticity point in 
the stiffness. The increase of stiffness is related to the 
incompressible nature of the insulating phases 4^ For the 
system we studied, the phase is of Luther-Emery type for 
small />2£ and a composite phase of Luther-Emery-liquid 
type in the wings and insulting type in the bulk for large 
p. We test for a system in a confinement of strength 
V2/t = 10"'^ and Ns = 200 lattice sites. The comparison 
of SDFT with DMRG shows that, in Fig. El the stiffness 
calculated from SDFT is quantitatively correct for weak 
interaction strength oi u ~ —1 and qualitatively correct 
for u — —2. The large deviation happens for the sys- 
tem of large number of particles. We find by DMRG, 
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TABLE 1: Ground-state energies per site (in units of t) with 
liarmonic trap Kxt/i = 4 x 10^^, and Kxt/t = 1 x 10^* for 
the system of = 100, Nf = 30 and iV^ = 300, Nt = 90 
respectively. The couple strength are u = —0.5, —1, —1.5, 
and —2. The SDFT results are compared to the DMRG ones. 
\AEo\ is the relative error between the SDFT and DMRG 
calculations. 
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Nr 
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30 
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different quantum phases 4^ In the top panel of Fig. [TOl 
the density-functional fidelity is plotted as a function of 
the coupling strength u with the distance between u and 
u' chosen as 0.002 and 0.004. We notice that the fidelity 
drops abruptly at the critical point u = —1.848, which re- 
fiects the fact that the two ground states in the different 
quantum phases have the maximum distance. In other 
words, when we fix the number of fcrmions and the am- 
plitude of the external potential by varying u, this critical 
value of the attractive interaction is obtained. At that 
point, the system is in the BI phase with the localized 
site occupancy of = 2 in the center of the trap. 

The fidelity susceptibility, the leading term of the fi- 
delity, is found to play a central role in the fidelity ap- 
proach in precisely characterizing the quantum phase 
transition, which is defined as, 

F{u,u + Su) = 1-^^XF, (18) 
where the fidelity susceptibility (FS) xf takes the form 

E l / dfiA'^ 
4n^ \duj ■ 

In the bottom panel of Fig. [TOl the FS is shown as a func- 
tion of the coupling strength u, which manifests itself as 
a marked peak at the same critical point as that in the 
fidelity function. The sharp peak in the FS implies prin- 
cipally that the density distribution evolves dramatically 
around the critical point, which is therefore taken as an- 
other important quantity to detect the quantum phase 
transition4^i^ The peak in the FS is free from the ar- 
bitrariness of the small parameter Su. To prove this, we 
calculated the FS with two different Su in Fig. [TOTb). 
Both cases have good agreements and support the con- 
clusion that the fidelity susceptibility does not depend on 
the value of Su, but the fidelity doesi^ 



the phase transition happens at A^f — 141 {Nf = 125) for 
u = —1 [u = —2), driven by the increasing number of 
the particles, while the results from SDFT slightly devi- 
ate from the DMRG ones with Af = 140 (Af = 124) for 
u = —1 (u = —2). We notice that the bulk compressibil- 
ity or the stiffness calculated here is not a real singular 
quantity since it includes the sizable contributions from 
the liquid wings of the density profilci^ 

Finally we study the phase transition driven by 
the attractive interaction strength through the density- 
functional fidelity, which is defined as^ 



F{u,u') = — tr^nj(it)ni(u') , (17) 
Nf 

here rij (u) is the density of particles in lattice site i with 
interaction strength u. If we fix the distance between u 
and w', Su = u — u' , the density- functional fidelity is ex- 
pected to show a drop around the critical point due to 
the maximum distance between the two ground states in 



V. SUMMARY 

In summary, in this paper we have studied the ground- 
state properties of two-component Fermi gases of attrac- 
tive interactions in a ID system, by performing SDFT 
and DMRG simulations. Three different systems have 
been tested, for fermions on a lattice line with periodic 
boundary condition, in a hard wall or in the presence 
of a harmonic confinement. For the system of periodic 
boundary condition or hard- wall, we showed the calcu- 
lations for the GSE and the energy gap in accordance 
with the exact Bethe-ansatz, CBCS, and DMRG results. 
We found that the ground-state energy is faithfully repro- 
duced by SDFT in the large range of interaction strength 
and fillings. For the system under the harmonic confine- 
ment, we calculated the pair-binding and four-particle 
binding energies to show that the system of attractive 
interactions induces effective attractions between parti- 
cles leading to pairs formation and effective repulsions 
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TABLE II: Pair- and four-particle binding energies £p and £4 per site (in units of t) are calculated with the energies from Table 
Ufor the system of A'f = 30 and TVs = 100 under the presence of a harmonic trap 14xt/t = 4 x 10"'^. The couple strength are 
u = —0.5, —1, —1.5, and —2. The SDFT results are compared to the DMRG ones. A£ is the relative error between the SDFT 
and DMRG calculations. 



u 


£p (SDFT) £, 


, (DMRG) A£j, 


£4 (SDFT) £4 (DMRG) Afi 


-0.5 


-0.000211 - 


-0.000207 1.9 


0.002054 0.002058 0.2 


-1 


-0.000559 - 


-0.000559 0.1 


0.001962 0.001952 0.5 


-1.5 


-0.001278 - 


-0.001219 4.8 


0.001935 0.001850 4.6 


-2 


-0.002881 - 


-0.002417 19.2 


0.001842 0.001754 5.0 



between pairs avoiding clusters formation. The system is 
either in the Luther-Emery type phase or in the compos- 
ite phase of Luthcr-Emcry-likc in the wings and insulat- 
ing like in the bulk. In order to characterize the quan- 
tum phase transition driven by the particle number or 
the attractive interaction strength, the thermodynamic 
stiffness, the density-functional fidelity, and the fidelity 
susceptibility are studied, respectively. We found that 
the thermodynamic stiffness is able to locate the quan- 
tum critical point driven by the particle number through 
the nonanalyticity point, and F (xf) is able to locate 
the quantum critical point driven by the attractive inter- 
action through the abrupt drop (the sharp peak) at the 
critical point. 
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